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Abstract 



The problem of confinement of fermions in 1+1 dimensions is approached 
with a linear potential in the Dirac equation by considering a mixing of 
Lorentz vector and scalar couplings. Analytical bound-states solutions are 
obtained when the scalar coupling is of sufficient intensity compared to the 
vector coupling. 



The Coulomb potential of a point electric charge in a 1+1 dimension, 
considered as the time component of a Lorentz vector, is linear and so it pro- 
vides a constant electric field always pointing to, or from, the point charge. 
This problem is related to the confinement of fermions in the Schwinger and 
in the massive Schwinger models 0-H and in the Thirring- Schwinger model 
||. It is frustrating that, due to the tunneling effect (Klein's paradox), there 
are no bound states for this kind of potential regardless of the strength of the 
potential ||]-||. The linear potential, considered as a Lorentz scalar, is also 
related to the quarkonium model in one-plus-one dimensions Recently 
it was incorrectly concluded that even in this case there is just one bound 
state 11. Later, the proper solutions for this last problem were found f|- 



TT|. However, it is well known from the quarkonium phenomenology in the 
real 3+1 dimensional world that the best fit for meson spectroscopy is found 
for a convenient mixture of vector and scalar potentials put by hand in the 
equations (see, e.g., ||12|). Therefore, the problem of confinement of fermions 
by a linear potential in 1+1 dimensions deserves a more general analyses. 
With this in mind, we approach in the present paper the Dirac equation 
in one-plus-one dimensions with a linear potential considering it as a more 
general mixing of Lorentz vector and scalar couplings. It is found that there 
are analytical bound-state solutions on condition that the scalar component 
of the potential is of sufficient strength compared to the vector component 
(|K| > l^tl)- As a by-product, the present approach also provides the oppor- 
tunity to find that there exist relativistic confining potentials providing no 
bound-state solutions in the nonrelativistic limit. Although we shall confine 
our discussion to the vector-scalar mixing, the inclusion of a pseudoscalar 
potential can be allowed. The hearth of the matter is a unitary transforma- 
tion similar to that one recently applied to the vector-pseudoscalar mixing in 
3+1 dimensions That unitary transformation leads to a Sturm-Liouville 
eigenvalue problem for the upper component of the Dirac spinor. 

Let us begin by presenting the Dirac equation in 1+1 dimensions. In 
the presence of a time- independent potential the 1+1 dimensional time- 
independent Dirac equation for a fermion of rest mass m reads 

Hty = (1) 



H = cap + (3mc 2 + V (2) 
where E is the energy of the fermion, c is the velocity of light and p is the 
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momentum operator, a and (3 are Hermitian square matrices satisfying the 
relations a 2 = j3 2 = 1, {a, f3} = 0. From the last two relations it steams that 
both a and (3 are traceless and have eigenvalues equal to —1, so that one 
can conclude that a and f3 are even-dimensional matrices. One can choose 
the 2x2 Pauli matrices satisfying the same algebra as a and (3, resulting in a 
2-component spinor The positive definite function \^\ 2 = satisfying 
a continuity equation, is interpreted as a probability position density and its 
norm is a constant of motion. This interpretation is completely satisfactory 



for single-particle states 0]. We use a = o\ and j3 = 03. For the potential 
matrix we consider 

V = lV t + (3V S + aV e + /3>fV p (3) 

where 1 stands for the 2x2 identity matrix and /?7 5 = <Ji- This is the most 
general combination of Lorentz structures for the potential matrix because 
there are only four linearly independent 2x2 matrices. The subscripts for the 
terms of potential denote their properties under a Lorentz transformation: t 
and e for the time and space components of the 2-vector potential, s and p 
for the scalar and pseudoscalar terms, respectively. 
Defining the spinor ip as 

ip = exp (I; A) * (4) 

where 



h 

A(x) = J X dx'^p- (5) 
the space component of the vector potential is gauged away 

p + — \ * = exp (^a) pip (6) 



c / 

so that the time-independent Dirac equation can be rewritten as follows: 

Hip = Eip (7) 

H = a lC p + a 2 V p + a 3 (mc 2 + V s ) + lV t (8) 

showing that the space component of a vector potential only contributes to 
change the spinors by a local phase factor. Introducing the unitary operator 
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U(9) = exp ( -i-j<?i 



(9) 



where 9 is a real quantity such that — 7r < 9 < tt, the transform of the 
Hamiltonian (Kl) takes the form 



H = UHU- 1 = a lC p + a 2 W p cos (0) - (mc 2 + V S ) sin (9) 



+CT3 



(mc 2 + V.) cos (9) - V p sin (9) 



(10) 



In terms of the upper and the lower components of the transform of the 
spinor ip under the action of the operator U: 



X 



and, moreover, choosing 



V t = V s cos (9) - ^sin(0) 
the Dirac equation decomposes into: 



(11) 



(12) 



-h 2 c 2 6" + i he 



V s ' sin (9) + cos (0)] + [(mc 2 + V s ) sin (9) + V p cos (9) 

(13) 

+2 \E + mc 2 cos (9)] [V s cos (9) - V p sin (9)] - [e 2 - m 2 c 4 cos 2 (9)] } = 



_ . -hecj)' + [(mc 2 + V s ) sin (9) + V p cos , , . 
+ mr cos (9) 

where the prime denotes differentiation with respect to x. In terms of the 
upper and the lower components the spinor is normalized as 

f + \ X \ 2 )dx = l (15) 
so that both <fi and x are square integrable functions. 
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Now, we shall restrict our discussion to linear potentials with V p = 0, 
namely 



a s \x\ 



(16) 



v t = v s cos{e) 

i.e., \V S \ > \Vt\. Thus, Eg . Cp~3|) becomes the Schrodinger-like equation 



-h z c 



2 „2 j.// 



+ {a 2 sin 2 (fl)a; 2 + 2a s £cos (0) + mc 



(17) 



+e(x)hca s sin(0) - (# 2 - m 2 c 4 ) } = 



where e(ar) = for i / 0. It is clear that Eq.flTTD allows two distinct 

classes of solutions depending on sin(0). 
For the class sin(^) = 0, we define 



C± = <*±\x\ + ^ 



(18) 



where the plus sign corresponds to V t = V s (008(6*) = 1) and the minus sign 
corresponds to V t = —V s (cos(6 ) ) = — 1). Furthermore, 



± 



2a s 
h 2 c* 



E ± mc 2 



1/3 



2a, 



{e =F 



mc 



so that Eq.([r?D turns into the Airy differential equation 



(19) 



(20) 



which has square integrable solutions expressed in terms of the Airy functions 
T5f : 0(C±) = A±Ai(£±), where A± is a normalization constant. The joining 



condition of and its derivative at x = leads to the quantization conditions 
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Ai(/3±) = for odd parity solutions 

(21) 

Ai'{fi±) = for even parity solutions 
These quantization conditions have solutions only for (3± < and some of 



them is listed are Table I [15]. One can see that /3± < always corresponds 
to \E\ > mc 2 . Substitution of the roots of Ai(/3±) and Ai\(3±) into (|l]J) allow 
us to obtain the possible energies as the solutions of a forth-degree algebraic 
equation: 

E A =p 2mc 2 E 3 ± 2m 3 c e E - [m 4 c 8 + (2Tica s f |/3±| 3 ] = (22) 

Instead of giving explicit solutions to this algebraic equation in terms of rad- 
icals, we satisfy ourselves verifying that the roots of Eq. (E2p always satisfy 
the requirement \E\ > mc 2 by using the Descartes' rule of signs (hence- 
forth DRS). The DRS states that an algebraic equation with real coefficients 
dk\ k + ... + a\X + cio = the difference between the number of changes of 
signs in the sequence a^, ai, do and the number of positive real roots is an 
even number or zero, with a root of multiplicity k counted as k roots and 
not counting the null coefficients (see, e.g., |I7||-||IB||). The verification of the 



existence of solutions for E > mc 2 is made simpler if we write E = mc 2 + 5. 
We get 



S 4 + 2mc 2 5 3 - (2hca s y 



5 4 + Qmc 2 5 6 + 12m 2 cV + 8mV5 - (2hca s y = 



(23) 



Observing the difference of signs among the coefficients of the leading coef- 
ficient and the lowest degree it becomes clear that there exist positive roots. 
In this particular case the DRS assures that there exists just one solution, 
since there is only one change of sign in the sequence of coefficients of (P3|). 
It is interesting to note that this result is true whatever the fermion masses 
and the coupling constant. The very same conclusion for E < mc 2 can be ob- 
tained by observing that the upper-sign solutions (V t = V s ) are mapped into 
the lower-sign solutions (V t = —V s ), and vice- versa, by the change E — > —E. 
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On the other hand, for the class sin(#) 7^ 0, corresponding to \V t \ < \V S 
we define 

£= — (\ x \ + x o) 

x 

where 



E cos (9) + mc 2 

x = 



a s sin (9) 



(24) 



/2a s sin(. yi 
^ = V he X ° 



Moreover, 



[E + mc 2 cos (9)} - m 2 c 4 cos (29) 
+ 2ftca s sin 3 (#) ( } 



such that 



d 2 m , e 



2 



(1/+ 1/2)^(0 a;>0 



f/( , + ^(0 = (26) 

whose solutions are the square integrable parabolic cylinder functions |15| : 
0(0 = BD V (£), for x > 0, and 0(0 = for a; < 0. B and C are 

normalization constants. Making use of the recurrence formulas 



D' v (z) - Z -D v (z) +D V+1 (z) = 



D' v (z) + ^D u (z)-vD v ^(z) = 



(27) 



and the matching conditions at x — 0, the quantization condition is 



D u+1 (£0) = iv^+LD, (£0) (28) 

where £0 is given by (f24|). By solving the quantization condition (|28D for 
z/ imposing that the solutions of ( p5| ) vanish for £ — > +00, one obtains the 
possible energy levels by inserting those allowed values of v in 
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E = —mc 2 cos(#) ± \Jm 2 c A cos {26) +2hca s {u + l)sin 3 (fl) 



The numerical computation of (|28| ) is substantially simpler when D v+ \ (£o) 



is written in terms of Z)[, (£o): 

"6 



^ (Co) 



9 =F 1 



(fo) (29) 



Because the normalization of the spinor is not important for the calculation 
of the spectrum, one can arbitrarily choose D v (£ ) = 1. By using a fourth- 
fifth order Runge-Kutta method |L6] an infinite sequence of allowed values of 
v are found corresponding to each sign in (|29"D . The lowest states are given 
in Table II for £o — 1- 

It is worthwhile to note that the first class of solutions (|Vs| = \Vt\) is 
independent of the sign of a s whereas the second class of solutions (|Vg| > 
|VJ|) depends on \a s sin(^)|. This observation permit us to conclude that even 
a "repulsive" potential can be a confining potential. This peculiar effect 
is because the scalar potential behaves like an x-dependent rest mass [pTf . 
Nevertheless, only potentials with a s > in the mixing a t = a s cos(6 l ) are 
confining potentials in the nonrelativistic approximation and the case V t = 
—V s reduces, for any sign of a s , to the case of a free fermion. It is well known 
that a confining potential in the nonrelativistic approach is not confining 
in the relativistic approach when it is considered as a Lorentz vector. It 
surprising that relativistic confining potentials may result in nonconfinement 
in the nonrelativistic approach. This last phenomenon is a consequence of the 
fact that scalar and vector potentials couples differently in the Dirac equation 
whereas there is no such distinction between scalar and vector potentials 
in the Schrodinger equation. Therefore the results exposed in this paper, 
beyond to present a generalization of previous results, might be of relevance 
to the quarkonium phenomenology in a four-dimensional space-time. 
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Table 1: The first roots of the Airy function and its first derivative 



M-l/fel) = o 


Ai(-\/3 ± \) = 


1.019 


2.338 


3.248 


4.088 


4.820 


5.521 



Table 2: The lowest solutions of Eq. (28) for £ = 1- 



z/ for the minus sign 


v for the plus sign 


1.580xl0" 4 
2.681 
5.038 


1.525 
3.915 
6.210 
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